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D | Abstract. We consider the problem of counting the number of 3-colourings of the 

edges (bonds) of the 4-8 lattice and the 3-12 lattice. These lattices are Archimedean 
with coordination number 3, and can be regarded as decorated versions of the square 
and honeycomb lattice, respectively. We solve these edge-colouring problems in the 
• infinite-lattice limit by mapping them to other models whose solution is known. The 

colouring problem on the 4-8 lattice is mapped to a completely packed loop model with 
loop fugacity n = 3 on the square lattice, which in turn can be mapped to a six-vertex 
I model. The colouring problem on the 3-12 lattice is mapped to the same problem on 

the honeycomb lattice. The 3-edge-colouring problems on the 4-8 and 3-12 lattices 
are equivalent to the 3-vertex-colouring problems (and thus to the zero-temperature 
Q \ 3-state antiferromagnetic Potts model) on the "square kagome" ("squagome") and 

"triangular kagome" lattices, respectively. 
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1. Introduction 

Counting problems involving strong local constraints arise in many different contexts. 
One important class of such problems deals with graph colourings, which is a topic 
of interest in several fields, including combinatorics and graph theory pQ, theoretical 
computer science, and statistical mechanics. In the latter field, connections have been 
identified between various colouring problems and other types of problems, such as Potts 
models [2], loop models [3111], and tiling [5] and folding [6] problems. Recently quantum 
generalizations of some colouring problems have also been studied [7]. 

If each vertex of a graph G is coloured with one out of q colours, such that no vertices 
connected by an edge have the same colour, the result is said to be a q- vertex-colouring 
ofG0 The number Pg{q) of such colourings for a graph with N vertices is a polynomial 
in q of degree N, called the chromatic polynomial of G, and the smallest integer q for 
which Pciq) > is the chromatic number xg [I]- The problem of determining Pg{q) 
for a given q is the q- vertex-colouring problem on G. For q > Xg, Pg(q) is equal to the 
zero-temperature partition function of the g-state antiferromagnetic (AF) Potts model 
on G. These AF Potts models are, for sufficiently large q, examples of nn-frustrated 
systems with an extensive ground state entropy: the ground state degeneracy increases 
exponentially with N, while the energy associated with each edge takes the lowest value 
possible. 

Similar considerations apply to colourings of the edges of a graph. If each edge of a 
graph G is coloured with one out of q colours, such that no edges meeting at a vertex have 
the same colour, the result is said to be a g-edge-colouring of G. The number of such 
colourings is the edge chromatic polynomial P^(g), and the edge chromatic number 
X*q is the smallest integer q such that P G e \q) > 0. The (/-edge-colouring problem on G 
is equivalent to the q- vertex-colouring problem on a related graph G', defined such that 
every edge of G hosts a vertex of G', and two vertices in G' are connected by an edge 
if their host edges in G meet at a vertex. Nevertheless, in spite of this equivalence, it 
may still be preferable to study this problem in its original edge-colouring formulation. 

An interesting class of lattices in two dimensions are the Archimedean lattices, 
which consist of uniform tilings of the plane with regular polygons such that all vertices 
have an identical environment. There are 11 such lattices, the most familiar ones being 
the square, honeycomb, triangular, and kagome lattices. The edge chromatic number for 
an Archimedean lattice is equal to its coordination number [8]. The simplest nontrivial 
g-edge-colouring problem for an Archimedean lattice is therefore obtained by taking q 
equal to its coordination number. This is the problem we consider in the following. 
Let £ be a lattice with iV vertices and coordination number q, with periodic boundary 
conditions in both directions, so the system is topologically a torus. The number of 
g-edge-colourings Zc scales as Wj? as iV — > oo. That is, Wc is defined as 



W c = lim Z 



l/N 




TV-s-oo 



Vertices and edges are often instead called sites and bonds in the physics literature. 
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Previously, Wc has been calculated exactly for the 3-edge-colouring problem on the 
honeycomb lattice [3] and for the 4-edge-colouring problem on the square lattice [9]. 
In these works the problem was formulated in terms of a transfer matrix that was 
diagonalized with a coordinate Bethe Ansatz [10] . 

In this paper we calculate Wc exactly for the 3-edge-colouring problem on the 
4-8 and 3-12 lattices (12] (see Fig. [T|), both of which are Archimedean lattices with 
coordination number 3. Another common feature of these lattices is that they can both 
be regarded as a "decorated" version of a "parent" lattice; for the 4-8 (3-12) lattice the 
parent lattice is the square (honeycomb) lattice. We show that these edge-colouring 
problems can be mapped to previously solved models on the "parent" lattices. For the 
4-8 lattice the mapping is to a completely packed loop model on the square lattice, which 
in turn can be mapped to a six-vertex model. For the 3-12 lattice the mapping (which 
is very simple in this case) is to the 3-edge-colouring problem on the honeycomb lattice. 
We also briefly discuss the alternative formulation of the 3-edge-colouring problems on 
the 4-8 and 3-12 lattices as 3- vertex-colouring problems (and thus zero-temperature 3- 
state AF Potts models) on the "square kagome" ("squagome") and "triangular kagome" 
lattices [13], and we comment on the 3-edge-colouring problem on the 4-6-12 lattice, 
another Archimedean lattice with coordination number 3. 



(b) 



Figure 1. (a) The 4-8 lattice, (b) The 3-12 lattice. 



2. The 3-edge-colouring problem on the 4-8 lattice 

2.1. Identifying the possible local 3-colouring configurations 

The 4-8 lattice is shown in Fig. [Ha). It can be regarded as a "decorated" square lattice, 
the decoration procedure consisting of replacing every vertex on the square lattice with 
a square rotated by 45 degrees. We will denote the edges of squares on the 4-8 lattice as 
"inner" edges and the edges connecting neighboring squares as "outer" edges. Consider 
now one of these squares together with the four outer edges touching it, i.e. a collection 
of eight edges in total. By considering all possible ways of 3-colouring the four outer 
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edges, and then for each such way considering 3-colourings of the four inner edges, we 
find that there are only 18 possible (i.e. allowed) 3-colourings of the eight edges. In 12 
of these colourings, each of the horizontal outer edges is "paired" with a vertical outer 
edge in the sense that edges in the same pair have the same colour and edges in different 
pairs have different colour. For each such colouring of the outer edges there is only one 
way to colour the inner edges. Examples of these colourings are shown in Fig. |2l top 
row, (a)-(d). In the remaining 6 colourings, all four outer edges have the same colour. 
For each colour of the outer edges there are then two possible ways to colour the inner 
edges. Examples of these colourings are shown in Fig. [2j top row, (e)-(f). (Note that 
the 12 colourings not shown in Fig. [2] can be generated by cyclically permuting the 3 
colours in the 6 colourings that are shown.) 




(a) (b) (c) (d) (e) (f) 

Figure 2. (in colour) Top row: 6 of the 18 possible local 3-colouring configurations 
on the 4-8 lattice. Bottom row: the corresponding local coloured CPL configurations 
on the square lattice. The remaining 12 configurations can be obtained by cyclically 
permuting the 3 colours in the configurations shown here. We have here ordered the 
colours cyclically as blue — > red — > green — > blue (see text for further details). 



2.2. Mapping to a completely packed loop (CPL) model on the square lattice 

The partition function of the completely packed loop (CPL) model with loop fugacity 
n on the square lattice is given by 

c 

Here the sum is over all configurations of completely packed loops, and M c is the 
number of loops in configuration c. By definition, a CPL configuration consists of 
non-intersecting loops such that every edge on the square lattice is covered by one loop 
segment [H]. This leads to two possible arrangements of loop segments meeting at a 
vertex, as shown in Fig. |3j 

The 3-edge-colouring problem on the 4-8 lattice can be mapped to a CPL model 
with loop fugacity n = 3 on the square lattice. The mapping is established by showing 
that each of the 18 possible local 3-colouring configurations can be mapped to a local 
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Figure 3. The two possible CPL configurations at a vertex on the square lattice. 

CPL configuration of the types shown in Fig. [3j but with loops that can take 3 possible 
colours (note that the number of such coloured local CPL configurations is indeed 
18 = 2 • 3 2 ). Let us first consider the 12 colourings with two different colours on 
the outer edges. By forgetting about the inner edges (as noted above, their colouring is 
uniquely determined in this case and thus represents redundant information) and simply 
connecting outer edges with the same colour, we get CPL arrangements in which the 
two meeting loop segments have different colours. This is illustrated in Fig. E](a)-(d). 
On the other hand, for the 6 colourings in which the outer edges have the same colour, 
there are two possible ways to colour the inner edges, which again can be mapped to 
the two ways that loop segments in the CPL model (now of the same colour, namely 
that of the outer edges) meet at a vertex. For example, one can make the rule that 
in this case two outer edges are to be connected if the inner edge between them in the 
3-colouring has a colour that follows cyclically after the colour of the outer edges. This 
is illustrated in Fig. [2] (e)-(f). Thus we arrive at a model of completely packed loops 
in which each loop can take 3 possible colours. It follows that ^4_8, the number of 
3-edge-colourings on the 4-8 lattice, equals the number of loop configurations in this 
coloured CPL model. This latter number is in turn equal to the partition function of a 
CPL model of un-coloured loops in which each loop is given a weight of 3, i.e., 

Z 4 _ 8 = Z(3). (2.2) 

We note that the CPL model (12.11) is in a non-critical phase for n = 3 (the model is 
critical for \n\ < 2) [14J. 

2.3. Mapping the CPL model to a six-vertex model. Calculation ofW^s 

The CPL model can be mapped to a six- vertex model [T5J [16]. To this end, we first 
rewrite Zin) as a sum over configurations of directed loops. Writing n = oj + a; -1 , a 
directed loop is assigned a weight u ±l where the two signs correspond to the two possible 
directions of the loop. Next one writes these directed loop weights as products of local 
weights associated with each vertex visited by the loop: the assigned vertex weight is 
Co> 1//4 or u; -1 / 4 depending on whether the directed loop turns right or left at the vertex. 
There is however a subtlety here [T5], [16] . While the mapping to local vertex weights 
gives the correct loop weight uj + u' 1 = n for all contractible loops, on a torus one can 
also have non-contractible loops (i.e. loops that wind around the torus), and these get 
an incorrect weight of 1 + 1 = 2 from this procedure. This corresponds to a modified 
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loop model with partition function 

Z{n) = ^n^2 M ^\ (2.3) 

c 

where A/" c c (A/" c nc ) is the number of contractible (non-contractible) loops in configuration 
c. We then have 

Z(n) = Z(n) ■ ^y^y (2.4) 

where the expectation value is taken with respect to the model defined by (12. 3p . Let 
the square lattice have M x (M y ) vertices in the x (y) direction, i.e. M = M x M y vertices 
in total. Without loss of generality we take M y < M x . The maximum number of 
non-contractible loops a loop configuration can possibly contain is then M x . Such a 
configuration would have no contractible loops, while the M x non-contractible loops, 
each of length 2M y , would wind around the torus in the y direction. This implies 
that (n/2) Mx is an upper bound for the expectation value in (12.41) . Furthermore, 1 is 
clearly a lower bound. These bounds imply that for the 4-8 lattice with N = 4M sites, 
W^s = liniAr^oo zV3, satisfies 



lim Z(3) 1/(4M) < W A _ 8 < lim 

M->oo M-5>oo 



2(3) 



2 \ M x 



1/(4M) 



lim Z(3) 1/(4M) (2.5) 

M->oo 



2 / 

(here limM^oo = lim^Af^oo)) from which it follows that 

W 4 - 8 = lim Z(3) 1/(4M) . (2.6) 

M->oo 

Zin) can be interpreted [TBI [TB"] as the partition function of a six- vertex model 
on the square lattice with periodic boundary conditions in both directions and model 
parameters given by 

a = 6 = l, c = oo 1 / 2 + oo- 1 / 2 = Vn~+2 A = - + 9 ~ ^ = --. (2.7) 

2ab 2 

Thus A = —3/2 for our colouring problem, so the six- vertex model is in the non-critical 
"antiferroelectric" phase that occurs for A < — 1. In this phase the free energy per 
vertex (site), /, is given by p2] 



+ J2 — ^7m- < 2 - 8 > 

< ' m man m A 



kuT 2 ^— ' m cosh(mA) ' 



m=l 

-A 



where the parameter A > is defined by A = — cosh A. Introducing x = e we find 
that 



e 



fc=i v 7 

For A = — 3/2, x = 2 — <p where ip = (1 + \/5)/2 is the golden ratio, which implies 
x 1/2 (l - x) = x. Thus W 4 _ 8 = ( e - //fesT ) 1/4 is given by 

^-" 1/4 n& = 1-24048... (2.10) 

fe=i v ' ' 
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3. The 3-edge-colouring problem on the 3-12 lattice 

The 3-12 lattice is shown in Fig. QJb). It can be regarded as a "decorated" honeycomb 
(he) lattice, the decoration procedure consisting of replacing every vertex on the 
honeycomb lattice with a triangle that points in opposite directions (up or down) on the 
two sublattices. Completely analogous to what we did for the 4-8 lattice, we denote the 
edges of the triangles as "inner" edges and the edges connecting neighboring triangles as 
"outer" edges, and then we consider the possible 3-colourings of the edges of a triangle 
and the three outer edges touching it. One easily sees that the only allowed 3-colourings 
are such that the three outer edges all have different colours, in which case the colouring 
of the inner edges is unique (namely the colour of an inner edge must be the same as that 
of the outer edge with which it does not share a vertex) [Fig. B^a)]. We can therefore 
forget about the inner edges and so a local 3-colouring on the 3-12 lattice maps to a local 
3-colouring on the original honeycomb lattice formed by connecting the outer edges as 
shown in Fig. 11(b) . Thus, as the number of 3-colourings is the same on the two lattices, 
while the number of vertices differ by a factor of 3, we get 

W 3 _i 2 = w£ /3 = (^] y/T(l/3) = 1.06522 . . . , (3.1) 



where we made use of the exact result for Wh c [3]. 






(a) (b) (c) 

Figure 4. (in colour) (a) On the 3-12 lattice, the three outer edges touched by the 
same triangle must all have different colours. Given the colours on the outer edges, 
the inner edges can only be coloured in one way. (b) Mapping to a 3-colouring on the 
honeycomb lattice, obtained from (a) by removing the inner edges and connecting the 
outer ones, (c) (See |Appcndix A3| ) A dimer covering on the 3-12 lattice that contains 
triangles without dimers (as shown here; edges occupied by dimers are red, unoccupied 
edges are black) can not be made into a 3-colouring, because the fact that the triangle 
edges form an odd-length loop makes it impossible to satisfy the colouring constraints 
when they are coloured with the remaining two colours. 



4. Final remarks 

As noted in Sec. HJ a g-edge-colouring problem on some lattice can be reformulated 
as a g-vertex-colouring problem on a different but related lattice. In particular, the 3- 
edge-colouring problem on the 4-8 (3-12) lattice is equivalent to the 3-vertex-colouring 
problem on the "square kagome" ("triangular kagome") lattice [TB]. The number of 
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3-vertex-colourings is in turn equal to the number of ground states in the 3-state AF 

Potts model on the same lattice. Defining s as the ground state entropy per vertex of 

the Potts model then gives (sk/tk = square/triangular kagome lattice) 

2 2 
s s k = g log W4-8, stk = g log W 3 -i2- (4.1) 

It is known that some edge-colouring problems are equivalent to packed loop 
(FPL) models (or variants thereof) on the same lattice. The 3-edge-colouring problem on 
the honeycomb lattice can be mapped to an FPL model {3], while the 4-edge-colouring 
problem on the square lattice can be mapped to a fully packed double loop (FPL 2 ) 
model [1]. The 3-edge-colouring problem on the 4-8 lattice can also be mapped to an 
FPL model [see Eq. flA.ip ]; the mapping is identical to that for the honeycomb lattice. 
In these mappings the loops have fugacity n = 2. In Ref. [18J it was shown that the 
n = 2 FPL model on the 4-8 lattice can be mapped to a particular case of a square- 
lattice 0(2) model that is further mappable to the n = 3 CPL model [19]. This route 
via an FPL and 0(2) model therefore constitutes an alternative way of mapping the 
3-edge-colouring problem on the 4-8 lattice to the n = 3 CPL model on the square 
lattice. The direct mapping given in Sec. 12.21 is however simpler. 

Another example of a colouring problem that can be mapped to a CPL model is 
the 3-edge-colouring problem on the 4-6-12 lattice, which is another Archimedean lattice 
with coordination number q = 3. It can be regarded as a "decorated" version of the 
kagome lattice, in which each vertex of the kagome lattice has been replaced by a square, 
in exactly the same way as when creating the 4-8 lattice from the square lattice (like the 
square lattice, the kagome lattice also has coordination number 4). Thus "locally" the 

3- colorings are the same as for the 4-8 lattice. Therefore it follows from the analysis in 
Sec. 12.11 and 12.21 that there are 18 possible local 3-edge-colouring configurations on the 

4- 6-12 lattice as well, and the 3-edge-colouring problem on this lattice can be mapped 
to a CPL model on the kagome lattice with loop fugacity n = 3. 
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Appendix A. On bounds for Wc 

Appendix A.l. The 4-8 lattice 

It is instructive and useful to check the exact result ( 12 . 1 0[) for W^s against some lower 
and upper bounds. We first consider a lower bound that can be derived from the CPL 
model on the square lattice. Since for positive n, Z{n) is a monotonically increasing 
function of n, one can get a lower bound for W^s by replacing Z 4 _ 8 = Z(3) with Z{1) 
in the definition (II. ip . Note that Z{1) is just the total number of CPL configurations on 
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the square lattice, which equals 2 M since there are two possible local CPL configurations 
(Fig. [3]) which can be selected independently at each vertex. Since M = N/4 this gives 
a lower bound of 2 1 / 4 rs 1.189 for W^s- 

An alternative and (as it will turn out) slightly better lower bound for W±-% 
can be established by regarding an arbitrary 3-edge-colouring on the 4-8 lattice as a 
configuration of fully packed loops on this lattice. To see this, imagine trying to make 
a 3-edge-colouring by first colouring all edges that have, say, colour A. These edges 
constitute a dimer covering of the 4-8 lattice since each vertex must be touched by an 
A-coloured edge. The remaining edges will form loops on the 4-8 lattice. At each vertex, 
two edges are covered by a loop while the third edge is "empty" (i.e. has colour A). 
As each vertex is visited by one loop, this is a configuration of fully packed loops [T4] . 
Since the 4-8 lattice is bipartite, all loops will cover an even number of edges, and so 
it is possible to colour neighboring edges along each loop in alternating colours as • • • 
B-C-B-C-B-C • • •, thus indeed giving a valid 3-colouring. As each loop can be coloured 
in two different ways (which are related by a shift of the colours by one edge along the 
loop), the number of 3-colourings is given by 



where the sum is over all dimer coverings of the 4-8 lattice and Md is the number of 
loops formed by the edges that are not covered by a dimer. Thus a lower bound for Z 4 _g 
is 34-8, the number of dimer coverings on the 4-8 lattice. It follows that a lower bound 
for VU 4 _ 8 is given by 2U 4 _ 8 = lim^^^, which equals [2D|I2T] Vl-457... « 1.207. 

The B-coloured edges also constitute a dimer covering, and so do the C-coloured 
edges. Thus a 3-edge-colouring can be regarded as consisting of 3 non-overlapping 
dimer coverings, one of each colour [3j. When two non-overlapping dimer coverings 
have been deposited, the third and last one is completely determined. Therefore an 
upper bound for Z 4 _ 8 is given by 3 4 _ 8 , the number of ways of depositing two dimer 
coverings completely independently [3]. Thus SU^.s = 1.457. . . is an upper bound for 



A slightly better upper bound can be obtained from 1 1A.1I) by replacing Md with 
its maximum over the set of all dimer coverings d of A-coloured edges. This maximum 
is obtained when the loops of B- and C-coloured edges are as short as possible (since 
the total length of all loops is the same for all dimer coverings). This occurs when 
each square is covered by such a loop (i.e. the A-coloured edges are then the "outer" 
edges connecting neighbouring squares). Thus Md = N/4 for this dimer covering, giving 
Z 4 -s < 3 4 - 8 • 2 W/4 . It follows that 2B 4 _ 8 * 21/4 = L435 • • • is an upper bound for W 4 _ 8 . 

Using the best lower and upper bounds established above, we have therefore found 
that 1.207... < IU 4 _ 8 < 1.435..., which is consistent with the exact result 1 12. 10p . 
We see that the exact result is quite close to the lower bound; the ratio between 
these numbers is ~ 1.027. This ratio is very similar to the corresponding one for the 
honeycomb lattice, which like the 4-8 lattice is bipartite with coordination number 
3: Using the results in Refs. [3 J and [20J, the ratio for the honeycomb lattice is 




(A.l) 



a 



W 4 . 



8 
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^1.460... /1.381... « 1.028. 
Appendix A. 2. The 3-12 lattice 

Defining 2U3_i2 = \ivsx^^ > . 00 ^)\_ vi where 33-12 is the number of dimer coverings on the 
3-12 lattice, we have [201 122] 2B 3 -i 2 = 2 1 / 6 « 1.122 which is greater than the exact 
result for W3-12 in ( 13. ip . Thus for the 3-12 lattice, the number of dimer coverings is not 
a lower bound for the number of 3-colourings, unlike the situation for the 4-8 lattice. 
This is because on the 3-12 lattice, there are dimer coverings that cannot be used as a 



starting point for generating 3-colourings (in the sense discussed in Appendix A.l ). This 
is related to the fact that the 3-12 lattice is not bipartite, which leads to the existence 
of odd-length loops, as illustrated in Fig. Htc) . 
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